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Abstract. We employ the Planck 2013 CMB temperature anisotropy and lensing data, and 
baryon acoustic oscillation (BAO) data to constrain a phenomenological rcCDM model, where 
dark matter and dark energy interact. We assume time-dependent equation of state param¬ 
eter for dark energy, and treat dark matter and dark energy as fluids whose energy-exchange 
rate is proportional to the dark-matter density. The CMB data alone leave a strong de¬ 
generacy between the interaction rate and the physical CDM density parameter today, Wc, 
allowing a large interaction rate |r| ~ Hq. However, as has been known for a while, the 
BAO data break this degeneracy. Moreover, we exploit the CMB lensing potential likeli¬ 
hood, which probes the matter perturbations at redshift z ~ 2 and is very sensitive to the 
growth of structure, and hence one of the tools for discerning between the ACDM model 
and its alternatives. However, we find that in the non-phantom models (rcde > ~1)) the 
constraints remain unchanged by the inclusion of the lensing data and consistent with zero 
interaction, —0.14 < T/Hq < 0.02 at 95% CL. On the contrary, in the phantom models 
(rcde < energy transfer from dark energy to dark matter is moderately favoured over 

the non-interacting model; —0.57 < T/Hq < —0.10 at 95% CL with CMB-I-BAO, while 
addition of the lensing data shifts this to —0.46 < T/Hq < —0.01. 
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1 Introduction 

The recently published Planck full-sky maps of the anisotropies of the cosmic microwave 
background (CMB) radiation are, in general, compatible with a simple 6-parameter cosmo¬ 
logical “standard” model [1-8] , the ACDM model, which was also in a good agreement with 
the earlier, less accurate, measurements. Apart from the 5% ordinary matter and radiation 
components, two exotic building blocks are needed. 68% of the energy density of the Universe 
comes from dark energy and 27% from dark matter. Dark energy causes the accelerating 
expansion of the Universe and dark matter is an invisible — weakly interacting or noninter¬ 
acting — component that is needed to explain the galaxy rotation curves [9], the structure 
formation, and the flatness of the spatial geometry of the Universe. Indeed, cosmological 
models with only ordinary matter do not provide a good fit to a broad variety of cosmologi¬ 
cal and astrophysical data. Thus, it can be regarded well established that some kind of “dark 
sector” is needed. 

Although the micro-physics of the dark sector is still largely unknown, many of its 
properties can be tested or constrained by comparing theoretical predictions to the CMB 
and large scale structure (LSS) observations. It has been proposed that the origin of both 
dark matter and dark energy could be some kind of scalar fields. From a particle physics 
point of view, it would be natural to assume that these fields interact with each other or 
with (dark) matter [10-43]. The dark matter could “decay” to dark energy and vice versa. 
However, in the cosmological standard model, dark matter and dark energy are assumed to 
feel only each others gravitational effects. In this paper we relax this hypothesis by allowing 
for a phenomenological energy transfer term, i.e., interaction between dark matter and dark 
energy. 

We treat the dark matter (c) and dark energy (de) as fluids that have equation of 
state parameters Wc = 0 and Wdeit) = Pde{t)/Pde{t) and whose energy-exchange rate is 
proportional to the dark matter density. While the total energy density obeys the usual 
continuity equation, 

Ptot 4 " 37^(1 -|- Wtot)ptot = 0 ) ( 1 - 1 ) 
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the model is characterized in the background by modified equations for the individual com¬ 
ponents, 


Pc + 3'Hpc = (a/aref)Qc, (1-2) 

Pde 37^(1 -|- Zn(Je)Pde — (®/®ref)Qde ; (1'3) 

where Qde = ~Qc; a prime indicates derivative with respect to conformal time r, a is the 
scale factor of the Universe, and Ti = a'/a is the conformal Hubble parameter. The scale 
factor at a reference time, Oref, is usually chosen to be today’s scale factor and normalized 
to unity, Oref = oq = 1, but for later convenience we keep it in the above equations. 

In the literature, there are three main categories for the choice of the phenomenological 
energy exchange term Qc'. 

(1) For example in [44-71] Qc assumed to be proportional to the Hubble rate, H, times 
either of the energy densities or their sum or some other combination of the energy 
densities. In [47] a phenomenological coupling of type Qc = constant xHpc is motivated 
by a quintessence model which induces a time dependent mass for the dark matter 
particles. The energy transfer could also change its direction during the history of 
the Universe, as studied in [72] by fitting to the background data an interaction of 
the form of a piecewise constant x Hp. The above-mentioned papers either focus only 
on the background evolution or, if perturbations are included, assume that H in the 
interaction term is spatially constant, i.e., assume it to be the average expansion rate. 
However, when deriving the perturbation equations in a consistent way, one should 
treat H as a local variable, i.e., include the term 6H [73], as also noticed in [74, 75]. 

(2) Another common choice for a phenomenological interaction is a constant times either 
of the energy densities or some combination of them (without including the Hubble 
parameter and hence an implicit time dependence), as done for example in [51, 56, 70, 
71, 76-91]. In other contexts, this type of interaction has been used to describe the 
decay of curvaton [92-95] or production of quintessence field condensation by a slow 
decay of superheavy dark matter [96]. 

(3) An interesting possibility is also an elastic scattering between the dark-sector “parti¬ 
cles”, which does not lead to the energy exchange in the background. These models 
leave much milder observational signatures and indeed the current upper limits on the 
interaction cross-section in these models are several orders of magnitude larger than 
the Thomson cross-section ax [97], if Wde > — 1- On the other hand, tight constraints 
(< 10“®iTr) are obtained if Wde < — 1 [98]. 

In this paper we will study a model where 

Qc = -Tpc, (1.4) 

with T = constant. A positive T describes an energy transfer or a decay of dark matter 
to dark energy {pc —>• pde in the background) and a negative T corresponds to an energy 
transfer from dark energy to dark matter (pde Pc in the background). From equations 
(1.2) and (1.3) we see that F has the same unit as H. Therefore, it is convenient to describe 
the interaction by a dimensionless constant T/Hq, which gives the interaction rate in units 
of today’s expansion rate. 
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The interaction (1.4) falls into the category (2) and we will follow closely [78, 81, 82], 
Our first aim is to update the constraints^ by using the Planck CMB temperature anisotropy 
data [2], including the full treatment of perturbations and calculation of the CMB angular 
power spectrum. As seen, e.g., in [82], the CMB data alone leave a strong degeneracy between 
the interaction rate and today’s physical CDM density parameter Wc, since the CMB probes 
the CDM density at the time of last scattering at redshift z* ~ 1100, and the interaction 
modifies the evolution between last scattering surface and today. Therefore, it is important 
to choose a complementary data set that constrains the matter (or dark energy) density 
near to today. In [82] we showed that either baryon acoustic oscillation (BAO) data or 
supernova data are almost “orthogonal” to the CMB data when constraining the interacting 
dark-sector model; see also [99] . To break the degeneracies we choose now three recent BAO 
measurements that probe the redshift range Zeff ~ 0.1 - 0.6: 6dF Galaxy Survey [100] at 
Zes = 0.106, Sloan Digital Sky Survey Data Release 7 (SDSS DR7) BAO [101] data point at 
Zes = 0.35 as reanalyzed by [102], and Baryon Oscillation Spectroscopic Survey Data Release 
9 (BOSS DR9) [103] at z^s = 0.57. These choices coincide with those made in the Planck 
2013 analysis in [3]. 

Among many other measurements than just the primary CMB temperature anisotropy, 
Planck made a more than 25a detection of weak gravitational lensing of CMB by large- 
scale structure [104] (which in the 2015 release improved to 4Qa [105]), superseding the 
previous ~ 5 ct detections of Atacama Cosmology Telescope (ACT) [106] and South Pole 
Telescope (SPT) [107]. The CMB lensing potential power spectrum, Cf^, probes the matter 
perturbations at redshift z ~ 2. It is very sensitive to the growth of structure and hence 
one of the tools for discerning between the ACDM model and its alternatives [108, 109], 
such as the interacting dark-sector models [75]. Thus, as a novelty and our second aim, we 
test the ability of the publicly available Planck 2013 likelihood to further constrain the 
interacting dark sector.^ 

The Planck temperature CMB data (without lensing likelihood) have already been used 
also in constraining the interacting models with scalar held dark energy [111], type (1) in¬ 
teractions Qc = [112-117], Qc = -h i2Pde) [118-120], Qc = H^pAePc/{Pde + Pc) 

[121], and type (2) interactions in [89] — but [89] neglects the perturbations and derives only 
background constraints, so it does not signihcantly overlap with our work. Yet a few other 
couplings are introduced in [122-126], with moderate Bayesian evidence for an interacting 
vacuum model [127], and steps toward more generic models are presented in [43, 128-131]. 

2 Background evolution 

As we will study interaction rates that are smaller than today’s expansion rate, \T\/Hq < 1, 
the effect on the evolution of pc is mild and until very recently pc oc a~^ to high accuracy. In 
contrary, the continuity equations (1.2) and (1.3) with our chosen interaction term (1.4) may 
lead to a negative dark-energy density, pde, at some point of the evolution of the Universe, 
see e.g. [78, 82]. As the nature and origin of dark energy is still a mystery, this behaviour 
could be acceptable. However, at the moment of zero crossing from positive values to the 

^In [82] only models with Wde > —1 were constrained. Now we stndy also phantom models {wde < —1) in 
separate Markov Chains Monte Carlo runs. 

^Although not explicitly clear, the Planck 2015 lensing likelihood (which has not been made publicly 
available as of 16th March 2015) may have been used in a different interacting dark-sector model in a Planck 
2015 paper [110], where a coupled dark energy model was constructed following [111]. 
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negative ones, the perturbation equations are singular [81]. Thus we will study only those 
models where pde is positive today and in the past. (It could cross zero in the future, but this 
is not a problem, since in the lack of observations from the future we do not know whether 
perturbations might start growing rapidly.) Here we briefly discuss the background evolution. 

We define a ratio of the energy densities, 


r = 


Pde 
Pc ’ 


( 2 . 1 ) 


and employ equations (1.2), (1.3) and (1.4) to find a differential equation 


r'= (—V+—T . (2.2) 

V ®ref / ®ref 

Had we studied a type (1) model with Qc = the resulting equation, r' = ^^(r — 

3wde)^ + "WT, would have been trivial to solve analytically. This mathematical simplicity 
seems to have been one of the reasons (although a questionable one) for the popularity of 
type (1) models in the literature. In our type (2) model, with Qc = —Tpc, the missing 
from the T terms in (2.2) makes full analytic solution impossible, and we need to solve the 
evolution of pc, Pdej and numerically. However, assuming radiation or matter domination, 
i.e., a known form for Ti, we can find analytic solutions. 

In the early universe, during radiation domination, Ti = and a/oref = t With 
the initial condition r(rref) = r^ei, we find 


r(r) 


= g(r/nef)(r2-T2^j)/2 


-) 

"Tef / 


-3lJ)de /..,.2 T \ -3«)de/2 

r„, + e(V'-)dV2) 

\ ^ T-ef / 


G 


+ 1, 2^ 


)-«( 


2^de T f 


r 

’ Uef 2 


(2.3) 


where G is the incomplete Gamma function, G{a,b) = x°’~^e~^dx. We choose Tref to be 

a later time than r, r^ef > r, and we assume that Pc(Tre{) and r^ef are positive. Now we 
figure out whether r, and thus pde, is positive at all r < Tref. The first term in (2.3) is clearly 
positive. The second term, which involves the difference of incomplete Gamma functions, 
can be written as 

(r/pef)(r/ 2 ) fz -) ( 2 . 4 ) 

V 2 TrefxJ ^ ^ 

The term in front of the integral is always positive. If T <0, the integral is from a more 
negative number to a less negative number, while the integrand itself is positive (and real) in 
the whole range. Therefore, with T < 0 we have r(r) > 0 (at least) for 0 < r < r^ef. If T > 0, 
then the integral is from a larger positive number to a smaller positive number, while the 
integrand itself is again positive. Since the integration is to a “wrong” direction, the result 
is negative. It is more negative the larger the T is and the further back in time we go, and 
in the limit r —>■ 0 the integral diverges, if —3wde/2 > 1. Thus with T > 0 and Wde < —2/3, 
the integral approaches — oo, when r —?• 0. This means that for T > 0 we can find a positive 
solution for r, in the whole range 0 < r < Tref only if rcde > —2/3 and jT] is small enough. 

During matter domination ?{ = 2t~^ and a/oref = (T/rref)^, where Tref is now some 
later time at matter domination than r. Now we find 

( \ —6wd( 

— ] 

'^ref J 

G (2wde + 1, ^4") - G (2n;de + 1, 

\ 'ref / \ 'ref / 


^ ref “1“ ^ 


(r/r/^f)GV3) / ^ 


3 r. 


ref 


-^^de 
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(2.5) 








Wde < -1 

-1 < Wde < -2/3 

-2/3 < Wde <1/3 

F > 0 

negative at early times 

negative at early times 

positive, if F small enough 

F < 0 

positive 

positive 

positive 


Table 1. The sign of pde between t = 0 and today, when we assume that today pde, Pc and H are 
positive. Only models with positive pde are studied in this paper. Here Wde refers to the dark energy 
equation of state parameter at radiation and (early) matter dominated eras. 


Similar arguments as in the radiation dominated case show that for F < 0, r is positive. 
If r > 0 and —2t(;de ^ Ij he., Wde < ~l/2, then the term, which involves the difference of 
incomplete Gamma functions, would approach —oo as r —)> 0. However, we should remember 
that the matter dominated solution does not hold at the earliest times. So, if F is suitably 
close to zero, r may stay positive even if —2/3 < Wde < “1/2- 

In table 1 we summarize the background evolution. This behaviour is confirmed by 
our full numerical solutions. In the following, we parametrize the variable equation of state 
parameter of dark energy as rcde = wo — Wai^ — a) [132, 133], and rewrite it 

Wde = woa + WeCi-- a), (2.6) 

where We = wq + Wa denotes the early time value of rcde and wq the late time value. 

When analysing the results, one should bear in mind that, in the phantom case, only 
negative values of F are possible; see the first column of table 1. Although both negative 
and positive values of F are possible in the non-phantom case, the theoretical prior set 
by demanding positive definiteness of the past background-energy densities leads to much 
larger parameter-space volume of models with a negative F. As the last column of table 1 
indicates, in the non-phantom models with positive F we must have Wde > —2/3 at early times 
and simultaneously |F| small enough, in order to have the required background behaviour. 
Moreover, the perturbations ’’blow-up” if We < —0.8 in the non-phantom models [81]. This 
adds another parameter-space volume effect to the analysis. The combination of these two 
effects is demonstrated in figure 1. 

3 Perturbation evolution 

Since the model is phenomenological, the perturbation equations do not follow directly from 
variation of a given action. We need to do some choices and keep in mind that these can 
not be completely arbitrary as the formulation should be covariant. A good example is type 
(1) models where H should not be treated as an average expansion rate, but instead 6H 
should be taken into account. For interaction (1.4) plausible perturbation equations were 
first derived in [78]. We will treat the perturbations in the same way as in [78] and [81], 
and assume a constant sound speed, = 1, for the dark energy. Here we only repeat the 
end result for the evolution of the dark energy and dark matter density contrast and their 
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Figure 1. The We volume effect in the non-phantom models {wde > ~1)- The points indicate 
samples from our Monte Carlo Markov Chains with the data combination CMB+BAO. The colour 
scale shows the value of wq for each sample. Whilst marginalizing (integrating) over the We direction, 
the models with a positive T (energy transfer from CDM to dark energy) receive much less weight 
than the ones with a negative F (energy transfer from dark energy to CDM). 


velocity potentials, 

^de + 3?^(1 — t(;de)5de + (1 + Wde) [^de + k^{B — E')'j + 9'H^ {1 — W1de)“^ T 


~ <5de + 3?^(1 — Wde)~j^ + ' 


aV 


Pc 


-3(1 -b Wde)'tp' = aT— 

Pde 

^de ~ “EHOde — J7~T ^'^de ~ ^ \ 

(1 + Wide) (1 + W^de) Pde 

+ 0c + k'^{B - E') - 31 / 1 ' = -aV(t ), 

e'^ + Wc -k‘^(j) = o, 


(0c - 20de) , 


A;2 

(3.1) 

(3.2) 

(3.3) 

(3.4) 


where 0, il), E, and B are scalar metric perturbations of spatially flat Friedmann-Robertson- 
Walker metric; 


ds^ = or 


I — (1 + 2(/))(ir^ -b 2diB drdx^ + (1 — 2'ijj)5ij + 2didjE dx^dx^'^ . (3.5) 


Along with the modified background evolution equations, we have implemented the per¬ 
turbation equations (3.1)-(3.4), into Code for Anisotropies in the Microwave Background 
(CAMB) in synchronous gauge where B = E = 0, and the velocity perturbations take the form 
0 = —k‘^v. In the last term on the first line of equation (3.1) we have, according to (2.6), 

wide = ('“^0 - We)a' = -Watt' = -WaOH. 


4 Data and methods 

To determine the posterior probability density function (pdf) of the parameters of our inter¬ 
acting dark-sector model, as well as the non-interacting model for a reference, we perform a 
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full Monte Carlo Markov Chain (MCMC) likelihood analysis, using our modified versions of 
CosmoMC [134, 135] and the CAMB [136] Boltzmann code. The data sets used in the analysis 
include the Planck 2013 temperature anisotropy [2] and lensing data [104], and three baryon 
acoustic oscillation (BAO) measurements [100, 102, 103], as described on page 3. In all the 
analysis we supplement Planck temperature data with the nine-year Wilkinson Microwave 
Anisotropy Probe (WMAP) polarization data [137], as done in [3], and label this combina¬ 
tion as “CMB”. For finding the best-fitting parameter values and the corresponding we 
use Bound Optimization BY Quadratic Approximation (BOBYQA) as provided in the cosmomc 
package, but with increased accuracy settings. 

In the MCMC scans and in the best-fit searches we assume spatially flat geometry 
of the Universe and vary the following cosmological parameters: the physical baryon and 
CDM densities today {iVh = and uoc = hP'VLc)-, the Hubble parameter, i.e., the expansion 
rate today {Hq in kms“^ Mpc“^), the optical depth due to reionization (r), the dark energy 
equation of state parameter today and at early times {wq and We), and the spectral index and 
logarithm of the amplitude of primordial scalar perturbations (ug and InA^). From these 
we can calculate various derived parameters, such as the dark energy density parameter 
today (Hde) and the current age of the Universe (Age). In addition we vary the nuisance 
parameters of the Planck likelihood code as per [3], and in the interacting dark-sector case 
also the interaction rate in units of today’s expansion rate, T/Hq, in the range (-1,1). 

5 Results 

We start with the non-phantom interacting model, which turns out not to improve the fit to 
the data over the non-interacting model. Then we proceed to the more interesting case of 
phantom models where the interaction can moderately improve the fit to the data, and the 
interaction can significantly affect the inferred values of standard cosmological parameters. 

5.1 Non-phantom model {w^e > ~1) 

In figure 2 we represent 68% and 95% confidence level (CL) contours in 2d slices of various 
standard cosmological parameters versus the interaction rate. The current CMB data alone 
leave almost as long degeneracy line in {ujc, P /Hq) plane as already found in [82] with the 
WMAP five-year [138] and the Arcminute Cosmology Bolometer Array Receiver (ACBAR) 
[139] data. This is natural, since CMB probes Wc*, i.e., the physical CDM density parameter 
at last scattering, not the density today. A positive F implies that the CDM density decreases 
faster than in the non-interacting case. So, in order to fit well to the acoustic peak structure 
of the CMB data, i.e., to have the same uJc* in the interacting model (with a positive F) as in 
the non-interacting model, we need a smaller ojc today. On the other hand, with a negative 
F we need a larger Wc today, in order to have a correct Wc*. 

However, we immediately notice also some differences to [82]. Most notably, the Planck 
data disfavour positive values of F over the negative ones more strongly than the WMAP 
data (see, e.g. the second panel of figure 2, the fourth panel of figure 3, and compare also 
figure 1 to a similar one in the appendix of [82]). A reason for this is that a negative F leads to 
a larger ojc, and hence to a smaller (derived parameter) Hde- This smaller Hde in turn causes 
much reduced integrated Sachs-Wolfe (ISW) effect compared to the non-interacting model 
or to the model with a positive F. At the lowest multipoles {i ~ 2-40), the Planck data have 
less temperature angular power (compared to the high multipoles) than the non-interacting 
ACDM model predicts. As the ISW effect adds power to the lowest multipoles the interacting 
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Figure 2. The non-phantom interacting model {wde > ~1)- 68% and 95% CL regions with 
CMB (gray), CMB+BAO (red), and CMB+BAO+lensing (blue) data. 
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Figure 3. Comparison of the non-phantom models (rcde > ~1)- Id marginalized pos¬ 
terior probability densities for the interacting dark-sector model (solid lines) and for the non¬ 
interacting “standard” rcCDM model (dashed lines). In the non-interacting model, with the CMB 
and CMB+BAO data, we allow the early-time dark energy equation of state parameter, Wg, to vary 
between -1 and +1/3, whereas with the CMB+BAO+lensing data we allow only the same range, 
Wg S (—0.8,+1/3), which is possible in the interacting model (see hgure 1), in order to provide a 
different perspective for comparing the models. 


model with a negative T can significantly improve the fit in this range, whereas a positive T 
leads to a worse fit. 
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CMB CMB+BAO CMB+BAO+lensing 


WMAP polarization (lowlike) 

2.16 

1.24 

0.00 

-0.03 

Planck \ow-l TT, 2 < ^ < 49 (commander) 

-1.24 

-1.32 

-0.10 

-0.16 

Planck high-.^ TT, 50 < I" < 2500 (CAMspec) 

-2.02 

0.67 

1.51 

0.05 

Three BAO measurements [100, 102, 103] 

- 

1.06 

0.49 

0.04 

Planck CMB lensing {Cf^) 

- 

- 

-0.49 

0.00 

S 

-1.10 

1.64 

1.41 

-0.11 


Table 2. Difference of the the best-fit between the interacting and non-interacting non- 
phantom {wde > —1) model. The first three columns indicate the difference to the non-interacting 
model with — 1 < We < 1/3, while the last column (highlighted by bold face) gives the difference to 
the non-interacting model with —0.8 < < 1/3. Negative values mean that the interacting model 

provides a better fit to the indicated data combination than the non-interacting one. 

As already mentioned, some non-CMB data are needed to break the (wc, F) degeneracy, 
and a powerful choice is the BAO data; see, e.g., figure 4 in [82]. The BAO data would 
lead to almost orthogonal constraints compared to the gray CMB alone curves in the top-left 
corner in the second panel of figure 2, forcing the combined CMB-I-BAO good-fit region near 
to the zero interaction. Adding the (Planck 2013) CMB lensing data does not modify the 
constraints. 

In figure 3 we show Id marginalized pdfs for all the primary parameters (except r and 
In As) and for one derived parameter flje) comparing the results obtained in the interacting 
model to the ones in the non-interacting model. This comparison is not straightforward, since 
the CMB data (and in particular CMB-I-BAO or CMB-|-BAO-|-lensing) would favour We ~ 
— 1, which is excluded in the interacting model (recall figure 1 which depicts the theoretical 
prior on We)- The consequence is that in the non-interacting model the most favoured values 
of We fall to the range —l<We< —0.8, which is forbidden in the interacting model. In order 
to demonstrate more directly what kind of effects the interaction has on the determination 
of the other parameters, we have made the MCMC run with CMB-|-BAO-|-lensing data in 
the non-interacting model with a prior —0.8 < We < 1/3. Hence, the solid blue (T ^ 0) 
and dashed blue (T = 0) lines indicate best the differences between the interacting and 
non-interacting model: in the interacting model a smaller Hq and fide and larger cOc are 
allowed (actually by any combination of data studied here). Table 4 summarizes the best¬ 
fitting values and 68% CL intervals of parameters, as well as, the shifts of the mean values 
of parameters compared to the non-interacting model in terms of the standard deviation (cr) 
of each parameter in the non-interacting model. For each parameter the first line represents 
the results in the interacting model and the second line in the non-interacting model. As 
was obvious from the figures, with CMB alone ujc shifts upward by more than 13 it and 
this is reflected in a shift downward by Hq {—2Aa) and Hde (—5.5 cj). Adding the BAO or 
BAO-blensing data, restricts these shifts to —1.6cr. The least sensitive parameters are r and 
InA^; these shift by less than 0.3 cj. 

Finally, we check whether the non-phantom interacting dark-sector model improves the 
fit to the data; see table 2. Once we break the (wc, F) degeneracy with BAO or BAO-|-lensing, 
the interacting model actually worsens the fit to the data compared to the non-interacting 
model, if we allow in the non-interacting model —l<We< 1/3. This is due to the fact 
that the combinations CMB-I-BAO or CMB-|-BAO-|-lensing strongly favour Wg < —0.8, as 
already mentioned. In the last column we allow for the non-interacting model only values 
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CMB 

CMB+BAO 

CMB+BAO+lensing 

WMAP polarization (lowlike) 

1.48 

0.88 

0.95 

Planck \ow-l TT, 2 < ^ < 49 (commander) 

-1.14 

-2.25 

-2.10 

Planck high-.^ TT, 50 < I" < 2500 (CAMspec) 

-1.05 

-2.37 

-0.41 

Three BAO measurements [100, 102, 103] 

- 

-0.14 

-0.13 

Planck CMB lensing {Cf^) 

- 

- 

-0.47 

S 

-0.71 

-3.87 

-2.17 


Table 3. Difference of the the best-fit between the interacting and non-interacting 
phantom (wde < —1) model. 


We > —0.8, and in this case we naturally find that the interacting model with its one extra 
parameter compared to the non-interacting model indeed improves the fit slightly, Ax^ = 
—0.11. The largest improvement, Ax^ = —0.16, comes from the Planck \ow-l TT data due 
to the decreased flje and thus a reduced ISW contribution. The overall conclusion is that 
the data do not favour the non-phantom interacting model. 

5.2 Phantom model (tCde < —1) 

In the phantom models the interacting and non-interacting cases are directly comparable, 
since all the standard parameters can have the same prior in both cases. In particular, we 
choose the uniform priors —3 < tco < — 1 and —3 < We < —1. Now the comparison of 
of the best-fitting interacting and non-interacting models shows a clear preference for 
the interacting model; see table 3. With CMB+BAO+lensing data the largest improvement 
comes again from the \ow-^ Planck TT data, where the ISW effect can be reduced in the 
interacting model. 

In the phantom model the 2d slices reveal many interesting things; see figure 4. Firstly, 
with CMB+BAO data, the non-interacting model is outside the 95% CL region, except in the 
(r, P) and {wq, F) figures, where F = 0 is at the border of this region. The marginalized Id pdf 
of F /Ho in the fourth panel of figure 5 thus shows a clear peak at as high an interaction rate as 
T/Hq = —0.4, about 3a away from zero. According to table 5 even the mean, F /Hq = —0.34, 
is almost 3cr away from zero. Therefore, with CMB+BAO data we might conclude that they 
moderately favour the interacting phantom model, giving an improvement of 3.9 for the best 
fit x^ with only one extra parameter. 

Adding the CMB lensing data to the constraint budget, restores the non-interacting 
model to or near to the border of the 68% CL region; see the blue curves in figure 4. We 
underline that there are two “separate” effects in the game: (a) An indirect improvement of 
the constraints for F due to the improved constraints of the equation of state parameters wq 
and We, which is visible also in the non-interacting model — compare the sequence of the 
gray, red, and blue dashed curves in figure 5 {wq and We panels). In the interacting model, 
adding BAO (red solid) to CMB alone (gray solid) does not change the pdf of We almost 
at all, but the addition of the lensing data (blue solid) leads to a solid constraint on We 
even. Moreover, figure 5 and table 5 reveal that (both in the non-interacting and interacting 
models) adding the lensing data shifts ojc to smaller values than favoured by CMB+BAO. 
Recalling the (wc, F) degeneracy, this causes a shift to a less negative interaction rate, (b) As 
explained in Introduction, the lensing data constrain the growth rate of perturbations (mainly 
around redshift z ~ 2). Since the interaction directly affects the evolution of perturbations, 
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Figure 4. The phantom interacting model {wde < —1) with CMB, BAO, and lensing data. 


- CMB - CMB+BAO - CMB+BAO+Lensing 

CMB(r = 0) CMB+BAO (r = 0) CMB+BAO+Lensing (r = 0) 
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Wo We Os Qde 


Figure 5. The phantom model (w+e < ~1) with CMB, BAO, and lensing data. 


the lensing data hence probe also directly the interaction rate. The combined effect of (a) 
and (b) is to drive T closer to zero. However, the Id marginalized pdf of the interaction rate 
still peaks at a clearly non-zero value, T/Hq = —0.19 (hgure 5), the best-htting value being 
—0.12 and the mean being —0.23, which is about l.Gcr away from zero. Nevertheless, this 
demonstrates how powerful the lensing data are in discerning between the standard ACDM 
(or in our case rcCDM) and the alternative dark-energy or dark-sector models. 

Unlike in the non-phantom model, in the phantom model allowing for an interaction 
affects dramatically the preferred values of some standard cosmological parameters even with 
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CMB+BAO (or CMB+BAO+lensing) data. Naturally, the most affected is ujc, which shifts 
20a (15cr) upward. Table 5 shows that this is reflected in fide by a 4(T shift downward; from 
the mean value of ffde = 0.72 to 0.65 (0.66). The values of other parameters stay largely 
unaffected, except Ug which moves 2cj upward. Despite of this, the scale-invariant primordial 
spectrum remains 3.3a (4.5it) disfavoured by the data in the interacting model, while it is 
7a disfavoured in the non-interacting tcCDM model. 

6 Discussion 

We presented a numerical analysis of a phenomenological interacting dark-sector model, 
where dark energy was treated as a fluid with a time-dependent equation of state parameter, 
tCde = woa + We{l — a). We considered a model with energy transfer from dark matter to 
dark energy of the form Qc = —7'pc, and assumed the relative energy transfer rate T to 
be a constant in time. We updated the previous analysis [82] of the non-phantom model 
Wde > —1, using the Planck 2013 data, and — for completeness — extended the analysis also 
to the phantom dark energy, rcde < — 1. For the first time, in the analysis of the interacting 
dark sector, we utilized the CMB leasing potential reconstruction data {Cf'^). 

In the non-phantom models the non-interacting tcCDM model (T = 0) lied close to 
the peak of the posterior pdf; —0.069 < T/Hq < 0.005 at 68% CL and —0.14 < T/Hq < 
0.02 at 95% CL, with CMB-|-BAO-|-lensing, and allowing for the interaction (and thus one 
extra parameter) did not improve the fit to the data. On the contrary, in the phantom 
models an interaction rate higher than 10% of the expansion rate, T/Hq = —0.16 (—0.12), 
led to an improvement of by 3.9 (2.2) with CMB-I-BAO (CMB-|-BAO-|-lensing) data. 
The non-interacting model lied outside of the 95% CL interval, —0.57 < T/Hq < —0.10 
(—0.46 < r /Hq < —0.01). Our results in the phantom model indicated that the lensing data 
have a great potential for discerning between interacting and non-interacting models, and 
considerably improve the constraints on T, as well as, on the dark energy equation of state 
parameter today and at early times, wq and We- 

We focused on the effects that can be calculated within linear perturbation theory (by 
using our modified version of CAMB), and thus we used only CMB and BAO data, and in 
particular demonstrated the discerning power of the CMB lensing data. Although the shape 
and amplitude of the matter power spectrum might provide some extra constraints, the 
use of it up to the smallest scales (largest wave numbers k) probed by the current large- 
scale structure observations would require detailed modeling of the non-linear effects that 
may differ from the non-interacting models, as shown in [140-144]. Yet there would remain 
ambiguity of the assumptions made on the bias function, i.e., how the dark matter traces the 
visible matter on the smallest observable scales. Going to even smaller scales would require, 
e.g., the use of Lyman-a forest data [145]. 

In the near future, the Dark Energy Survey (DES) [146] and later the Euclid [147] 
data will offer an unprecedented discerning power. Dark-sector interaction can increase the 
probability of very massive clusters occurring at high redshift [148] and existence of structures 
at very large scales [149]. Voids in the interacting models are emptier and contain less 
neutral hydrogen than in the ACDM model, and this effect might alleviate tensions between 
simulations and observations in voids [145]. Further, weak lensing has been shown to offer 
valuable constraints in the future [75, 150]. In this paper, we mentioned several times the 
modified ISW effect in the interacting dark-sector models. More direct ISW data (other than 
an increased \ow-^ temperature angular power which is impossible to distinguish from other 
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effects) is provided by the cross-correlation of the CMB and large-scale structure maps; see, 
e.g, [151-154]. As already mentioned, an interaction also directly modifies the growth rate of 
matter perturbations [155], which can make it difficult to distinguish from modified gravity 
[156]. In the background level, the effects of interaction may be absorbed into a modified 
effective equation of state parameter r(;de,eff; see [78]. Hence an interacting model might seem 
like a non-interacting model with a different equation of state parameter. This degeneracy is 
studied in [157]. Extensive summaries of the constraining or distinguishing power of various 
existing and future low-redshift (and) large-scale surveys are provided by [158, 159]. 

In this paper we adopted a phenomenological “top-down” approach, i.e., introduced an 
energy(-momentum) transfer term to the background (and perturbation) continuity equa¬ 
tions and studied whether such an extension of the rcCDM model was favoured by the data. 
Another possibility would have been a “bottom-up” approach, i.e., identifying fundamental 
theories where an interaction arises and then studying a particular form of the interaction 
predicted by the theory. Recently, an interesting scenario, which leads naturally to an in¬ 
teracting dark sector, has been put forward [160]. The so called disformal coupling (see also 

[130] ) arises from Dirac-Born-Infeld actions in Type II string theories, when matter resides 
on a moving hidden sector D-brane. As a bonus, the coupling in this scenario can alleviate 
the coincidence problem of dark energy. Since the current data do not exclude interactions in 
the dark sector and indeed can accommodate quite large interaction rates, this theory offers 
interesting prospects for the future work. 

Recently, progress in phenomenological or effective description has also been made by 
generalizing the parametrized post-Friedmann framework (PPF) [161] to the case of inter¬ 
acting dark energy [116, 162]. PPF can be used to describe very generic types of interaction 

[131] , at the cost of up to 12 free functions or extra parameters. However, when applied to 
specific models the number of extra degrees of freedom reduces dramatically. For example, 
in the case of our model to just one, as shown in [131]. 
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CMB 


CMB±BAO 

CMB±BAO±lensing 

Parameter 

Best fit 

68% limits 

Best fit 

68% limits 

Best fit 

68% limits 

II 

[-H 

3 

0) 

0.022169 

0.022096 

0.02212 ± 0.00029 (±0.3cr) 
0.02204 ± 0.00029 

0.022200 

0.022121 

0.02235 ± 0.00027 (±0.6cr) 
0.02219 ± 0.00026 

0.022298 

0.022275 

0.02238 ± 0.00026 (±0.3ct) 
0.02229 ± 0.00025 

UJc 

(r = 

0) 

0.1716 

0.11971 

0.155l)i;;;^(; (±13.2a) 

0.1199 ±0.0027 

0.11788 

0.11912 

0.1204+°:°°^^ (±1.4a) 

0.1177 ±0.0019 

0.11735 

0.11709 

0.12081°;°°^ (±2.5a) 

0.1165 ±0.0017 

II 

0) 

51.80 

67.06 

53.2i^;^ (-2.4a) 

66.73 

68.01 

64.9+)'i (-l.Ocr) 

66.2++° 

67.16 

67.20 

64.9+)l (-0.8cr) 

65.7l+®o 

T 

(r = 

0) 

0.0903 

0.0853 

0.089iroi^ 

0.0990 

0.0892 

0.097 ±0.014 (±0.3cr) 

0 093+° °^° 

u.uyo_Q g^4 

0.0924 

0.0932 

0.0991°°)® (±0.2a) 
0.0961°;°)® 

T/Ho 

(r = 

0) 

-0.475 

0 

-0.3410;)^ 

-0.0007 

0 

U.UdO_o.oi8 

-0.0024 

0 

u.udy_Q 020 

Wo 

(r = 

0) 

-0.835 

-1.000 

< -0.595 (±0.7cr) 

< -0.758 

-1.0000 

-1.0000 

< -0.946 (-0.5cr) 

< -0.911 

-0.9990 

-0.9998 

< -0.948 (-O.lcr) 

< -0.941 

We 

(r = 

0) 

-0.419 

-0.901 

-0.532+")'f (±1.6cr) 

< -0.737 

-0.7784 

-0.983 

-0.7041°;°°° (±2.4a) 

< -0.880 

-0.7987 

-0.7990 

-0.7061°;°°° (±0.7a) 

< -0.728 

ris 

(r = 

0) 

0.9620 

0.9630 

0.9611 ±0.0074 (±0.0ct) 

0.9608 ± 0.0074 

0.9670 

0.9634 

0.9703 ± 0.0066 (±0.8cr) 
0.9656 ± 0.0062 

0.9680 

0.9674 

0.9701 ± 0.0062 (±0.4cr) 

0.9676 ± 0.0060 

ln(10i°As) 

(r = o) 

3.0902 

3.0821 

3.090l(!;;;f7 (+0-1^^) 

3.089 ±0.025 

3.1045 

3.0879 

3.093 ±0.027 (±0.1cr) 
3.0901°;°®® 

3.0897 

3.0906 

3.098 ± 0.025 (±0.2ct) 
3.0941°;°®° 

^de 

(r = 

0) 

0.278 

0.6846 

0.36i^i" (-5.5a) 

0 626+° °®® 

U.OZU_o.o30 

0.6855 

0.6946 

0-66lll°^| (-1.4a) 

0.6801°;°)^ 

0.6904 

0.6914 

0.6591°;°)® (-1.6a) 

0.6781°;°)) 

Age/Gyr 

(r = o) 

14.415 

13.810 

14.33l°;^i (±3.3a) 

13.948;)!;)’f 

13.8243 

13.7785 

13.880 ±0.049 (±1.4ct) 
13.816l°;°g 

13.808 

13.8093 

13.8821°;°!) (±0.8a) 

13.8481°;°)^ 


Table 4. The non-phantom model {wde > ~1)- Best-fitting parameter values, mean and 68% CL intervals for selected cosmological parameters 
(the primary MCMC parameters in bold, and derived parameters in non-bold face). For each parameter the first line is for the interacting dark-sector 
model and the second line for the non-interacting model (F = 0). In parenthesis, the shift of the mean value compared to the non-interacting model 
is indicated in units of the standard deviation of the corresponding parameter of the non-interacting model. 






CMB 


CMB±BAO 

CMB±BAO±lensing 

Parameter 

Best fit 

68 % limits 

Best fit 

68 % limits 

Best fit 

68 % limits 

UJb 


0.022304 

0.02245 ± 0.00034 (-bl.lcr) 

0.022261 

0.02247 ± 0.00035 (±2.2cr) 

0.022332 


(r = 

0 ) 

0.022211 

0.02213 ± 0.00029 

0.022035 

0.02192 ± 0.00026 

0.022184 

0.02203 ± 0.00025 



0.1368 

0.162 ±0.023 (±16.0cr) 

0.1392 

“0l67™r7±20±+ 

0.1314 

0.1466+^-;)'(^" (±15.0(t) 

(r = 

0 ) 

0.11911 

0.1195 ±0.0026 

0.12063 

0.1223 ± 0.0022 

0.11864 

0.1201 ±0.0018 

Ho 


83.6 

81+i;; (-i.ia) 

69.30 

(±0.6cr) 

69.05 

70.8+^i (-0.4a) 

(r = 

0 ) 

98.6 

> 85.4 

69.42 

72.2l^-® 

68.88 

71.6tli 

T 


0.0892 

0.089ini4 (-O.Oct) 

0.0870 

0.089^u-u(^ (±0.3a) 

0.0860 

0.081 ±0.012 (±0.0 ct) 

(r = 

0 ) 

0.0928 

0.089iroll 

0.0766 

0.084 ±0.012 

0.0860 

0.080t°:°“ 

T/Ho 


-0.120 


-0.161 


-0.121 

— 0 9Q9+» J-4 

U.40Z_q ggy 

(r = 

0 ) 

0 

— 

0 

— 

0 

— 

Wo 


-1.53 

whole prior range (-3,-1) 

-1.100 

> -1.81 (-3.9cr) 

-1.129 

> -1.34 (-1.9ct) 

(r = 

0 ) 

-1.939 

-1.52 ±0.28 

-1.035 

> -1.19 

-1.021 

> -1.14 

We 


-2.18 

whole prior range (-3,-1) 

-1.50 

whole prior range (-3,-1) 

-1.043 

> -1.40 (-I.Ict) 

(r = 

0 ) 

-1.57 

whole prior range (-3,-1) 

-1.148 

> -1.35 

-1.007 

> -1.21 

ris 


0.9659 

0.9683 ± 0.0088 (±1.0cr) 

0.9640 

0.9691 ± 0.0094 (±2.2cr) 

0.9677 

+^9679™ri+1^6+ 

(r = 

0 ) 

0.9648 

0.9611 ±0.0073 

0.9594 

0.9550 ± 0.0064 

0.9640 

0.9581 ±0.0059 

ln(10^°As) 

3.0901 

3.0821^;^"27 (-0-2fT) 

3.0836 

3.081+^;^^^ (-0.2a) 

3.0755 

3.061 ± 0.023 (-0.5ct) 

(r = 

0 ) 

3.0959 

O.U50_o.o27 

3.0659 

3.085 ±0.024 

3.0799 

3.070j:°:°^° 

^de 


0.773 

0.695+;j;i|o (-3.7a) 

0.664 

0-647+";;j^« (-4.4a) 

0.6775 

o-eoiiU:"!? (-4.oa) 

(r = 

0 ) 

0.8548 

A Qi 7+0.043 

U.5i/_o.oi6 

0.7040 

A 7 Aq + 0.015 

U. <ZO_Q Q 20 

0.7032 

A 7 AA+O.OI 3 

U. /ZZ_0_017 

Age/Gyr 

13.580 

13.72+^;(^ (±2.7a) 

13.770 

13.792 ± 0.053 (±1.8ct) 

13.790 

13.789 ±0.054 (±1.6+ 

(r = 

0 ) 

13.423 

13.485;)!;)’^^ 

13.7474 

13.716 ±0.043 

13.7587 

13.720 ± 0.043 


Table 5. The phantom model (wde < ~1)- Best-fitting parameter values, mean and 68% CL intervals for selected cosmological parameters (the 
primary MCMC parameters in bold, and derived parameters in non-bold face). For each parameter the first line is for the interacting dark-sector 
model and the second line for the non-interacting model (F = 0). In parenthesis, the shift of the mean value compared to the non-interacting model 
is indicated in units of the standard deviation of the corresponding parameter of the non-interacting model. 
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